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Abstract 

We develop a detailed analysis of electron transport in normal diffusive con- 
ductors in the presence of proximity induced superconducting correlation. We 
calculated the linear conductance of the system and the profile of the electric 
field. A rich structure of temperature dependences was found and explained. 
Our results are consistent with recent experimental findings. 
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Recent progress in nanolithographic technology revived the interest to both experimental 
and theoretical investigation of electron transport in various mesoscopic proximity systems 
consisting of superconducting and normal metallic layers. In such systems the Cooper pair 
wave function of a superconductor penetrates into a normal metal at a distance which 
increases with decreasing temperature [jl| . At sufficiently low temperatures this distance be- 
comes large and the whole normal metal may acquire superconducting properties. Although 
this phenomenon has been already understood more than thirty years ago and intensively 
investigated during past decades, recently novel physical features of metallic proximity sys- 
tems have been discovered and studied theoretically (see P-P^ and further references 
therein). 

In this paper we study the influence of the proximity effect on transport properties 
of a diffusive conductor in the limit of relatively low temperatures and voltages. We will 
assume that this conductor is brought in a direct contact to a superconducting reservoir 
which serves as an effective injector of Cooper pairs into a normal metal. It is reasonable 
to assume that in this case the system conductance can only increase in comparison to 
its normal state value. As the penetration length of Cooper pairs into a normal metal 
becomes larger with decreasing temperature T one might expect the system conductance to 
increase with 1/T monotoneously. In what follows we shall present a rigorous calculation 
of a temperature dependent linear conductance of the the system and demonstrate that 
its actual temperature dependence turns out to be more complicated. In particular we 
will show that if the system contains no tunnel barriers there are two different physical 
regimes which determine the system conductance in different temperature intervals. At 
relatively high temperatures the superconducting correlation of electrons in the normal metal 



survives at the typical distance of order ~ y^/^) where V = VFUmp/^ is the diffusion 
coefficient. As T is lowered proximity induced superconductivity expands into the normal 
metal and, consequently, the "normally conducting" part of the system effectively shrinks in 
size. This effect results in increasing of the conductance of a normal metal. At sufficiently 
low temperature the length ^„ becomes of order of the size of the normal layer, i.e. the 
superconducting correlation survives everywhere in the normal metal and the effective gap 
in the quasiparticle spectrum develops there [|TB|. Physically this means that at sufficiently 
low energies and temperatures there exist no uncorrelated electrons in the normal metal 
with a superconducting reservoir attached to it. The characterisic energy which plays the 
role of such a gap is Sg ~ min(A, D/L'^). At energies lower than this value electron transport 
is only due to such correlated electrons (or Cooper pairs injected into the N-metall) which 
"density of states" increases with increasing e due to the gap effect. As the temperature is 
increased from zero higher and higher values of e contribute to the current and the system 
conductance increases with T at very low temperatures. Eventually with increasing T a 
crossover to a high temperature behavior of the system conductance takes place. Note 
that similar behavior of the normal metal conductance in the presence of proximity-induced 
superconductivity has been also found in a recent numerical analysis by Nazarov and Stoof 
^ . It is also interesting to point out that at T = the system conductance coincides with 



that of a normal metal with no proximity effects [|T^. 

We will also demonstrate that in the presence of tunnel barriers the temperature depen- 
dence of the conductance is essentially different. Whereas the contribution of the diffusive 
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region is nonmonotoneus in T, the conductance of the interfaces decreases monotonously 
with temperature. At the same time the proximity-induced superconducting gap developes 
in N [|I^. As a result, with lowering the barrier transparency the crossover takes place to 



the effective conductance decreasing monotonously with T, characteristic for two serial NIS' 
tunnel junctions (S' is now the diffusive normal conductor with proximity-induced gap). We 
note that both types of behavior, namely nonmonotonously enhanced and monotonously 
decreasing with T conductance have been observed in the experiment 

Let us consider a a quasi-one-dimensional normal conductor of a length 2L with a su- 
perconducting strip of a thickness 2ds attached to a normal metal on the top of it and two 
normal reservoirs attached to its edges (see fig.l). The length L is assumed to be much 
larger than the elastic mean free path limp but much shorter than the inelastic one. This 
geometrical realization has a direct relation to that investigated in the experiments [^Q. 
The two big normal reservoirs N' are assumed to be in thermodynamic equilibrium at the 
potentials V and respectively. In contrast to the case of ballistic constriction [113 the 



potential drop within the system is distributed between the interfaces and the conductor 
itself. The general approach to calculate the conductance of such kind of structures was 
developed in In what follows we shall apply this method to analyse the temperature 

dependence of the NS proximity structure of fig.l. 

The electron transport through the metallic system can be described by the equations 



V 



for a matrix of quasiclassical Green functionsG in the contact p!9|,p0 




(1) 



A A A 

where G^, and G^ are respectively the impurity-averaged advanced, retarded and 
Keldysh Green functions. These functions are in turn matrices in the Nambu space: 

A ^ A A AA/\/\A 

G'^=a, (?^+ tay -(G^)*and G'^=G'^f - fG^ 

A ^ A A 

. Here the distribution function /= /;+ a^/j, where // = tanh(e/2T) and ft describes 
deviation from nonequilibrium. Taking advantage of the normalization condition for the 
normal and the anomalous Green functions (g^)'^ — = 1 it is convenient to parametrize 
= cosh 6', = sinh6', where 6 = 6i + 162 is a complex function. Deep in the bulk 
superconductor it is equal to 9s = 1/2 In [(A + £:)/(A — e)] — in/2 for e < A and = 
(1/2) In [(e + A)/{e — A)] /2 for e > A (here and below we omit the indices R(A)). 

V 

The electrical current / and the electrostatic potential are expressed through G as 



de Sp a^G G 



(2) 



poo roo 

x)= deTig^ix.e) = / de ft{x,e)ve{x), (3) 
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where v is the the density of states, //^(x) = '^{gf{x)) and 5* is the crossection area of 
conductor. 

Being expressed in terms of the function 9{e,x) the equations p!9|j20|] for the Green 
functions and the distribution function for the N-metall take a particularly simple form 

Vdle + 2ie sinh ^ = (4) 



d. 



X 



Vicosh' 9,)dJt =0, (5) 



X is the coordinate along the N-conductor. Here we neglected the processes of inelastic 
relaxation and put the pair potential in the normal metal equal to zero A^y = assuming 
the absence of electron-electron interaction in this metal. 

Before we come to a detailed solution of the problem let us point out that the conclusion 
about the anomalous behavior of the system conductance can be reached already from the 
form of eq. @. Indeed it is quite clear from @ that the effective diffusion coefficient 
increases in the N-regions with proximity- induced superconductivity and, therefore, the 
electric field is partially expelled from these regions. This energy dependent field modulation 
is controlled by the solution for 6{e,x) and is directly related to the physical origin of the 
anomalous temperature dependence of the system conductance discussed below. 

The equations and should be supplemented by the boundary conditions at the 
interfaces of the normal metall N. Assuming that the anomalous Green function of big 
normal reservoirs N' is equal to zero from ||22|J^ we obtain 

iBdxO = sinh 6 , . 

jBCOshe.dJt = cosh9,{ft - ftix = 0,2L)), ^""^ 

where 'Jb = Rb/ Pn the interface resistance parameter, Rh is the resistance of the interface 
between the N-conductor and the N'-reservoirs, Ptv is the resistivity of the N- metal. In 
general we should also fix the boundary condition at the interface between the N-metal and 
the superconductor. For the case of a perfect transparency of this interface (which is only 
considered here) and for typical thickness of the normal layer wn ^ y/S 
C,N Cooper pairs easily penetrate into it due to the proximity effect and the Green functions 
of the N-metal at relatively low energies for d < x < d + 2ds are equal to those of a bulk 
superconductor 6 = Og. This intuitively obvious result can be proven rigorously (see e.g. 
|T5| and references therein). Thus the region of a normal metal situated directly under the 



superconductor for our purposes can be also considered as a piece of a superconductor S' 
and the solution of (H), (^ needs to be found only for < x < c? (without loss of generality 
we will stick to a symmetric configuration). 

Proceeding along the same lines as it has been done in ref. we arrive at the final 
expression for the current 

1 

I = — de [tanh((£ + eV)/2T) - tanh((£ + eV)/2T)] D{e), (7) 
2R Jo 

where D[e) defines the effective transparency of the system 

L + 'jB 



Die) 



cosh 92 (x=0,e) cos 61 (x=0,e) 



+ /q dxseoh 6i{x,e) 
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Let us consider the case of a sufficiently long normal conductor d ^ V/A. Then at low 
temperatures T -C A the interesting energy interval is restricted to e <^ A. For such values 
of e the contribution of the S"-part of the normal conductor shows no structure and can be 
easily taken into account with the aid of obvious relations 

I dx?,ec\i'9i{x,e) = / (ix sech^6'i (x, e) + (issech^6's i (8) 
Jo Jo 

and sech^6's 1 = ^1 — . Due to this reason we will discuss only the properties of the 
iV-part {0 < X < d). For the sake of completeness we will also demonstrate the effect of 
finite dg in the end of our calculation. 

For the differential conductance of the A^-part < x < d normalized to its normal 
( "non-proximity" ) value in the zero bias limit eq. (|^) yields 

-— =— / deD{e)sech\e/2T). (9) 



\dV J 2T Jo 

Analogously the normalized zero-bias electrostatic potential distribution reads 



^o[x) 



\im^ = ^["^ deD(e)u J x)sech\e/2T) [ dx' sech\9Ax')) (10) 
v^o V 2Td Jo Jx 



The analysis of the problem can be significantly simplified in the case of perfectly trans- 
parent interfaces (7^ = 0). In this case the boundary conditions are 

^(0) = (11) 
Oid) = Os (12) 

for the contact to the normal and the superconducting reservoir respectively. The effective 
transparency of the N-part then reads 

D{e)=(^j^dxsedi^{ei{x))^ . (13) 

As it was already pointed out for relatively long normal conductors and at low T only the 
energies e <C A give an important contribution to the conductance. In this case the typical 
energy scale is defined by the Thouless energy = "D/f^ <^ A. For these energies we can 
set 6s = iTx/2. Let us ffist put T = 0. Then the thermal distribution factor /5sech^(/3e) 
reduces to a delta function and we have 

(T^(T = 0) = D(0), (14) 
0o(a;) = (D(0)z/o(x)/c/) dx'sech^{ei{x')), (15) 



i.e. we only need the solution of (P with boundary conditions (|TT]) at £ = 0, which is 6* = i^x 
and therefore i'o{x) = cos(a;7r/2), where x = x/d. This means 



5 



ajv(T = 0) = 1. 



(16) 



This result coincides with that obtained first by Artemenko, Volkov and Zaitsev [|T^] and 
demonstrates that - in contrast to what one might expect intuitively - at T = and very 
low voltages the diffusive conductor does not "feel" the proximity-induced superconductivity 
and the conductance is exactly equal to its normal state value. The profile of the electric 
field penetrating into the normal metal at T = also can be found easily. We have 

Eo{x) := d^M^) = - 1) sin(x7r/2) - cos(x7r/2), (17) 

i.e. in the low temperature limit the electrical field distribution in the structure turns out 
to be essentially nonmonotoneous. In the case T -C we can calculate 6 perturbatively. 
From Vdf9 = — g)£ sinh 6'/(§) and (|lT]) we get 



6 = ^ — [x — sin(a;7r/2)] + i—x. 
£d 2 



Keeping only leading order terms in — , we get 



rp2 

CTTV = 1 + A — , (18) 

where A = ^ (| — C(2) ~ 0.049 is a universal constant. This means, that for low 
temperatures aN(T) grows quadratically on the scale of Ed and approaches the crossover 
towards the high temperature regime discussed below. 

We can also calculate the electrical field distribution in this approximation, however, 
the result is rather cumbersome and shows no structure which could not be seen from the 
numerical results presented below. 

In this limit T ^ (where we still have T ^ A), the contribution of the low energy 
components to the thermally weighted integrals for (Tn{T) and (poiT) is neglible and we only 
have to take into account the solutions of (^) for energies e 3> e^. It is well known (see e.g. 
|§]), that for this energy range the solution of (H) together with (|Tl|) reads 



tanh(^(x)/4) = tanh (^^) e'^^^"^) (19) 

where k = —2ie/V. By using obvious substitutions and multiple-argument relations for 
hyperbolic functions, we arrive at the following identity: 

/ dxsech{9^{x)) = il-x)~AJ^ / dy qiy) / (l + qiy)f (20) 

Jx y e Jo 

where q{y) = 4(3 + 2^/2)e~'^y sin^ y / {e-^y + 3 + 2^2)1 

For calculating D{e) we can, as the integrand becomes exponentially small for y > 
^{k) ^ 1, take the upper bound to infinity, such that it becomes a universal constant. 
From there we can calculate the conductivity in this limit 
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MT) = 1 + B^, (21) 



where again B = 0.42 is a universal constant. The whole integral in the rhs of (|20D is very 
small compared to one, so near the normal reservoir only the 1 — x-term has to be kept. As 
in this region u is, apart from exponentially small terms, constant, 0o is linear in x, so there 
Eq is constant. 

These results have a simple physical interpretation. Superconductivity penetrates into 
the normal part up to = \f^, whereas the rest stays normal, so the total voltage drops 
over a reduced distance d — ^n- Thus the resistance of the structure is reduced according to 
the Ohm law. In terms of the conductividy, this means 

= 1 + B'^ (22) 

which is equivalent to (pT]). 

For temperatures comparable to Ed the problem can only be handled numerically. The 
results show excellent agreement to the analytical calculations in the corresponding limits. 

The numerical results (see fig 2) confirm, that for -C A the universal scaling with ^ 
is excellently fulfilled, as even the peak in the conductivity is systematically a.t T ^ Se^ and 
the peak height of about 9% seems to be universal. This result also agrees with the recent 
numerical analysis |2l[] . This peak becomes smaller if we take into account the effect of finite 
ds keeping d fixed (fig. 3) The qualitative features, however, remain the same. 

For the field distribution at T <C we can observe the described nonmonotoneous 
behavior (see fig. 4). From T ^ Ed also a minimum in the distribution shows up. For T ^ Ed 
the electrical field is almost constant for x < 1 — ^at/c?, drops rapidly for a; > 1 — ^^Id and 
shows some interesting structure around x = 1 — ^^Id. 

Perhaps the most interesting effect found here is the nonmonotoneous dependence of the 
system conductance with temperature accompanied by a nontrivial shape of the electric field 
penetrating in the N-layer. As T is lowered proximity induced superconductivity in this layer 
becomes stronger and stronger being more and more efficient in expelling the field from the 
region adjacent to a superconductor. At very low temperatures this results in a substantial 
suppression of this field in a big part of the sample. Inevitably (the total potential drop is 
fixed!) this effect leads to an excess value of the electric field far from the NS boundary as 
it is shown in fig. 4. As it was already pointed out we can understand the effect of increase 
of the system conductance with increasing T at very low temperatures as a gap effect. At 
T Ed the conductance reaches its maximum and starts decreasing due to weakening of the 
proximity with further increase of T. 

Finally we would like to point out that the presence of tunnel barriers at NN' interfaces 
entirely changes the conductance of the system. If one lowers the barrier transparency the 
crossover takes place to the behavior demonstrating monotonously decreasing conductance 
with T (fig.5), which is typical for two serial NIS' tunnel junctions. Note that both types of 
behavior, namely nonmonotoneous and monotonously decreasing with T conductance have 
been observed in the experiment [Q. Further details will be presented elsewhere . 

We would like to thank G. Schon, C. Bruder and W. Belzig for useful discussions. This 
work was supported by the Deutsche Forschungsgemeinschaft within the Sonderforschungs- 
bereich 195. 
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FIGURES 
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FIG. 1. Physical setup. N' designates normal reservoirs, S superconductor, S' is effectly superconduct- 
ing due to the proximity to S, N is the normal part affected by the proximity to S' 
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Conductivity on the N part at e^«A 
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FIG. 2. Temperature dependence of the conductivity 
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FIG. 3. Influence of finite ds 
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FIG. 4. Field distribution for various temperatures 
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Conductivity in the presence of tunneling barriers, y=10 
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FIG. 5. Destruction of the structure of a by the presence of strong tunnehng barriers 
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